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Abstract 

Denote by A n the set of square (0, 1) matrices of order n. The set A n , n < 8, is 
partitioned into row/column permutation equivalence classes enabling derivation of 
various facts by simple counting. For example, the number of regular (0, 1) matrices 
of order 8 is 10160459763342013440. Let V n , S n denote the set of absolute determi- 
nant values and Smith normal forms of matrices from A n . Denote by a n the smallest 
integer not in T> n . The sets T>g and Sg are obtained; especially, ag = 103. The lower 
bounds for a n , 10 < n < 19, (exceeding the known lower bound a n > 2/ n _i, where 
f n is nth Fibonacci number) are obtained. Row/permutation equivalence classes of 
A n correspond to bipartite graphs with n black and n white vertices, and so the 
other applications of the classification are possible. 

Key words: (0, 1) matrices, Smith normal form, permutation equivalence, 
determinant range, classification 
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1 Introduction 



Let A n denote the set of square (0, 1) matrices of order n. Hadamard maximum 
determinant problem is: find the maximum determinant among the matrices 
in A n - In this paper we consider a slightly more general problem: determine 
the set V n = {\detA\ | A e A n }. 

It is known [1] that determinants of (0, 1) matrices of order n are related to 
determinants of ±1 matrices of order n + 1. If A is a (0, l)-matrix of order n, 
let B = ^(A) be a ±l-matrix of order n + 1 obtained from A by replacing 
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its by —1, bordering with a row — l's on the top, and a column of l's on 
the right. Clearly, \I/ is a one-to-one correspondence. By adding row 1 of B to 
each of the other rows of B, we see that det B = 2 n det A. 



By the Hadamard inequality | det B\ < y (n + l) n+1 , and therefore for all A e 

An | det A | < ^ n \J {n + l) n+1 . The equality is attained if B is an Hadamard 
matrix, i.e. if BB T = (n + l)I n+ i, where T denotes transposition, and I n is the 
unit matrix of order n; for n > 2 this implies n — Ak — 1. For upper bounds 
for determinants of A G A n see for example [2]. 

Let d n denote the largest element in V n , and let a n be the smallest integer not 
in V n . Craigen [3] shows that the set V n is the interval {1, 2, . . . , d n } for n < 6, 
but not for n = 7, because a 8 = 41 < d 8 = 56; he suggests that a 9 = 103. 

Some interesting sequences, related to (0, 1) matrices are found in [4]: A003432 
(the sequence d n ), A013588 (the sequence a n ), A051752 (c n , the number of 
matrices in A n with the determinant d n ) and A055165 (m n , the number of 
regular matrices in A n ). A few first members of these sequences are given in 
the following table. The values of a 9 , c 8 , Cg and m 8 seem to be new. 
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In this paper, which is a continuation of [5], the matrices in A n , n < 8, are 
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partitioned into row/column permutation equivalence classes, enabling the 
classification by ADV, and more precisely — by SNF (see section 2). Let S n 
denote the set of SNF's of matrices in A n . In section 3 the sets V 9 and £9 are 
determined. In section 4 the lower bounds for a n , 10 < n < 19 are obtained; 
c n , n < 9, are obtained in section 5. 

We introduce now some notation. If A = [a^] and B = [6^] are matrices of 
the same dimension m x n, we say that A < B if A is lexicographically less 
than B, i.e. if for some pair of indices (i,j) the first i — 1 rows of A and B 
are equal, the first j — 1 elements in the ith row of A and B are equal, and 
a-ij < hj. For example, 



1 




1 




< 




1 




1 1 



The smallest matrix in a set A is the representative of A. 

Denote by Pjj the permutation matrix obtained from /„ by exchanging the 
ith and jth row. 

The matrices A, B E A n are equivalent [7], A ~ B, if B is be obtained 
from A by a sequence of elementary row/column operations of the follow- 
ing types: exchange of two rows/columns, multiplication of a row/column by 
— 1, and addition/subtraction of a row/column to/from another row/column. 
Let SNF(A) denote the SNF of A. It is known that A ~ B is equivalent to 
SNF (A) = SNF(_B) (in [7] this statement is proved for polynomial matrices). 

The SNF diag(<ii, (I2, ■ ■ ■ , d n ) is written simply as a vector (di, d,2, ■ ■ ■ , d n ). If 
diagonal elements of SNF are repeated, we use the shortened exponential 
notation. For example, (1 3 ,2,0) is short (1,1,1,2,0). If s E S n , then we also 
say that the SNF-class s is the set {A E A | SNF(A) = s}. 

Let J n denote the square matrix of order n with all elements equal to one. 



2 Classification of (0, 1) matrices of order 8 or less 

The set V n could be obtained by computing determinants of all A E A n . A 
better approach is to group matrices with the same determinant, and then 
to compute the determinant of only one matrix in each group. It is useful 
to classify A n into subsets with constant absolute determinant value(ADV), 
or into even smaller subsets with constant SNF. We now review some such 
partitions of A n . 
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Let n r denote the group of row permutations of matrices from A n . Permuta- 
tions from Il r preserve ADV. 

The representative of the matrix A orbit is obtained from A by sorting its rows 
into a nondecreasing sequence. Rows of A correspond to binary numbers less 
than N = 2 n . Therefore, the number of orbits of Il r in A n is equal to (^"j -1 ) , 
i.e. the number of nondecreasing sequences of length n from {0, 1, . . . , TV — 1}. 
Let II denote the group of row and column permutations; II also preserves 
ADV. The group II induces an equivalence relation 7r over A n . We say that 
matrices A and B are permutationally equivalent, A ~ 7r B, if they are in the 
same orbit of II. Let A n denote the representative of the matrix A equivalence 
class (7r-class; we say shorter that A n is a ^-representative of A). 

Example 1 The it -representative of 

1 1 
1 1 
1 

is the matrix 

1 
Oil, 

1 1 

the smallest of all 36 permutationally equivalent matrices. 

Let A* denote the set of 7r-representatives in An. In [8] it is shown that the 
number of 7r-classes in A n is given by: 

\Al\= £ E C(i)C(i)exp 2 £ i r j s 2^ s \ (1) 

ii+2i2+...+ni n =n ji+2j2+---+nj n =n r,s=l 

where the summation is over all vectors i = (ii,i 2 , ■ ■ ■ ,i n ), j = (ji, J2, • • • ,jn), 
and 

C(i) =n\/{VH 1 \...n ln i n \) 

is the number n-permutations with i r cycles of length r, r = 1, 2, . . . , n; (r, s) 
denotes GCD of integers r, s. The values \A*\ are listed in Table 1; they are 
easily computed for quite a large n using, for example, UBASIC [9]. It is seen 
that p n is close to 2™ /(n\) 2 for n < 15. An effective algorithm to generate the 
representative A n of a given matrix A (section 2.3) simplifies the classification 
of matrices, because it enables to deal with the small subset A* of A n . 
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Table 1 



The number of permutationally nonequivalent matrices in A n , n < 15. 
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0.65875 


33642660 
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0.77266 


14685630688 
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0.85533 


21467043671008 


10 


0.91045 


105735224248507784 


1 1 

± ± 


94^6^ 


1 7643^69309^780761 4996 


12 


0.96754 


100455994644460412263071692 


13 


0.98088 


19674097197480928600253198363072 


14 


0.98886 


13363679231028322645152300040033513414 


15 


0.99358 


31735555932041230032311939400670284689732948 



2.1 Matrix extension 



In order to classify matrices in A n by ADV values, one has to select carefully 
the order by which determinants are computed. It is natural to start from 
matrices of order n — 1, and then to extend them by one row and one column 
of ones and zeros in each possible way. For an arbitrary B e A n -i, let bord(5) 
denote the subset of A n , containing matrices with the upper left minor equal 
to B. We say that the matrices in bord(-B) are obtained by extending B; if 
A G bord(fi), then A is an extension of B. 

The calculation of determinants of all matrices in bord(5) is an easy task. If 
A G bord(.B), then A is of the form 



x b 
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where x = [xi x 2 ... z„-i] and y = [ Vl y 2 ... y n -i] T . Then [1] 

n— 1 n— 1 

det A = 6 det -B — X! X! ^ (^) 

i=i i=i 

where is the cofactor of B, corresponding to a^. 
Obviously, 

A = {A | (5, y, 6) G A-i x {0, l}"- 1 x {0, l}"- 1 x {0, 1}}. 

If we precompute cofactors Bij, then determinant of each matrix from bord(5) 
is computed by only one addition: for the fixed x, the column y might traverse 
the set of possible values via a Gray code (so that in the sequence of y's each 
two subsequent vectors differ in exactly one position). 

Williamson [1] noted that it is enough to let B cross the set of 7r-representa- 
tives in A n -\. Let bord 7r (5) denote the set of ^-representatives of matrices in 
bord(5). 

Lemma 2 If B ~ vr B' then bord 7r ( J B) = bord^-B')- 



PROOF. Let A G bord n (B). If the row/column permutations, transforming 
B into B', are applied to the first n — 1 rows/columns of A, then the ma- 
trix with the upper left minor equal to B' is obtained. Therefore, the matrix 
permutationally equivalent to A is obtained by extending B', meaning that A 
is permutationally equivalent to a matrix from bord(-B'), i.e. A G bord n (B') . 
Analogously, bord 7r (i? / ) C bord„-(.B), and so bord 7r (5') = bord^-B). □ 



Not only determinants, but also SNF's of matrices in bord(5) can be ef- 
ficiently computed. The preprocessing step is to compute D = SNF(5) = 
diag(c?i, d 2 , ■ ■ ■ , d n ), and the matrices P, Q, such that PBQ = D, | det P\ = 
| detQ| = 1. In order to determine SNF(A) for an arbitrary A G bord(£?) of 
the form (2), we use the identity 

P 
1 

Denote xQ = [a\ a 2 ... a n ], Py = [c± c 2 ... c n ] T . Suppose d\ = d 2 = 
••• = dk = 1, for some k, 1 < k < n. Transforming the matrix from the 
righthand side of (4) by subtracting the row i multiplied by q from the row 
n, 1 < i < k, and then subtracting the column i multiplied by q from the 
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D Py 




x b 




1 




xQ b 



(4) 
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column n, 1 < % < k, we derive that A is equivalent to 
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Hence, SNF(A) determination is reduced to determination of SNF of a matrix 
of order n — k. The special cases when k > n — 1 are extremely simple, and 
they are not rare at all, because the corresponding SNF-classes are among the 
largest ones (at least for n < 9). More generally, one can reduce Oj, q modulo 
di, 1 < % < rank B. 

2.2 ^-extension 

Following Williamson [1], the approach based on extending 7r-representatives 
only, can be further improved. 

For an arbitrary A G A n let A' = X t A denote the matrix with the ith row 
equal to the ith row of A, and with the row j ^ i equal to the coordinatewise 
modulo two sum of jth and ith row of A. Equivalent ly, A' = RAS, where 
R is the matrix obtained from I n by subtracting ith. row from the others, 
and then by multiplying ith. row by —1; S is the matrix obtained from I n by 
changing sign of columns corresponding to ones in the ith row of A. A third 
equivalent definition of Xj [1] can be stated as follows: in the ±1 matrix 
B = ^(A) of order n + 1, the rows 1 and (i + 1) are exchanged, then the first 
row is "normalized" to all ones by changing signs of appropriate columns. By 
applying the matrix A' is obtained. Therefore, application of X t to A 
corresponds to a special row permutation in ^(A) (followed by scaling). It is 
natural to denote the identity transformation by X , X A = A. 

The transformation also preserves ADV. The composition of arbitrary two 
transformations X iy Xj is equivalent to only one: 

PijiXiA^tilKiJKntti^j 
A, if 1 < % = j < n 



Xi(XjA) = { 
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Let $ r denote the set of (n+ 1)! transforms of the form PXi, < i < n, where 
P is an arbitrary permutation matrix. 

Theorem 3 The set $ r is a transformation group of A n . 



PROOF. We have 



X t PA = PX Pi A, 



where pi is the index of the row of A, which is moved to the position i after the 
left multiplication by P. Let Pi and Pi be the two permutation matrices and 
let pj be the position to which Pi moves the row j after the left multiplication. 
Then 

_ { P2P1, Pj = * 

PiXjPxXi — P 2 P\X p Xi — < 

(/yyV< A'/<, />, / / 

If Pi = P is an arbitrary permutation matrix, 1 < i < n, P 2 = P~\ and 
Pj = i, then 

(px,)- 1 = p- l Xj. □ 



Clearly, each orbit of $ r contains at most n + 1 orbits of IT r . 

The corresponding transformation AXj over the columns of A (coordinatewise 
addition modulo two of the column i to all other columns) is defined by AXj = 
(XjA T ) T . Let $ c denote the group generated by column permutations and 
column transformations (-)Xi. 

Let $ be the group generated by the elements of groups $ r and $ c ; it also 
preserves ADV and its size is (n + l)! 2 . Matrices A and A' are said to be 
0-equivalent, A ~0 A', if they belong to the same orbit of $. Equivalently, 
A ~0 A' if and only if there exist row and column permutations P, Q, and 
row and column transformations Xi, Xj, such that A = PXiA'XjQ. For an 
arbitrary A e A n let A^ denote the ^-representative of A] 0-class of A is the 
orbit of $ containing A. 

Let bord^(P) denote the set of ^-representatives of matrices in bord(P). 
Williams [1] noted that $ and II have similar properties: in order to obtain the 
set A% of all ^-representatives in A n , it is enough to extend ^-representatives 
in An-i. 

Lemma 4 If B B' , then bord0(P) = bord^P/). 



PROOF. If B and B' are 0-equivalent, then there exist g G $, transforming 
B into B'. Suppose A e bord^(P). Then there exists a matrix A' e bord(P), 
A' ~0 A. By applying g to upper left minor of A, the matrix A" ~^ A', A" e 
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bord(-B') is obtained. Therefore, A ~0 A", and A G bord(-B'). Because A is a 
0-representative, we obtain A G bord^(-B'), implying bord^-B) C bord^(S'). 
Analogously, bord^-B') C bord^(-B), and hence bord^-B) = bord^(S'). □ 

2.3 Effective determination of it -representatives 

The classification of matrices in A n by extending matrices from »4^_i must 
be accompanied by an effective procedure to determine A n and A^ for an 
arbitrary A G A n . 

The matrix is the smallest among the family of at most n\ matrices obtained 
by sorting rows of all the column permutations of A. Search is performed more 
efficiently by a branch-and-bound algorithm. If we know the first % rows of A v 
(i.e. the row and column permutations P, Q such that the first % rows of PAQ 
are minimal), then the next row of A^ is a smallest column permutation (only 
permutations preserving the first i rows of PAQ are considered) of some of 
the remaining rows of PAQ. 

Algorithm 1 Branch-and-bound algorithm to determine A n given A G A n . 
Input : A G An 

Output : A n ; the permutation matrices P, Q, such that PAQ = A n ; 
count - the number of pairs (P,Q), such that PAQ = A n ; 

P(0) 4_ j • 0(0) T . A ^_ j . 

i <- 0; 
count <— 0; 
Optimize (i); 



{Continuation of the search for A w starting from the row i of p(* _1 ) AQ^ _1 \} 
{i.e. when the first i — 1 rows are already chosen and permuted} 
Optimize(i) 

Generate the minimal set of boundaries £W = (s^ = 0, Si \ . . . , sjj.'? = n) 
between adjacent columns of p( l_1 ) AQ^ % ~^ , such that the (i — l)-prefixes 
of columns from + 1 to are mutually equal, 1 < j < hi; 

for j = i to n do 

Vjl <- ErU,_ 1+ l (^"^^"^V 1 - Z - fci ' ^ nnm&er °/ 1 ' S } 

{in positions from + 1 to m £/ie jt/i row of P^'^AQ^ -1 ^} 

<— £/ie o/ indices of largest vectors Vj = (vji, Vj2, ■ ■ ■ , fjfcj, i < i < 7i; 
for a// j G do { £/ie candidates for the ith row of A n } 
p(i) ^_ p^.pG-i); { exchange the rows i and j} 
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compute Q^> from 1 \ so that all 1 's in the part of the row i from 

s[ t } 1 + 1 to sf ^ are moved to the right end of the part, 1 < I < ki; 

{hence preserving the first i — 1 rows of PW AQ^} 
compare the ith row of P^AQ^ to the ith row of A n : 
if the ith row of P^AQ® is less then 

copy the first i rows from P^AQ^ into A w ; 

fill with ones the rest of A n ; 

if i = n then P <— pW; Q <— ; count <— 1; else Optimize^ + 1); 
else if the ith row of P^AQ^ greater then 

continue; {bound step: try the next row index from L^} 
else 

if i — n then count <— count + 1; else Optimize^ + 1); 

Example 5 Algorithm 1, applied to the matrix from Example 1, gives the 
same n -representative as obtained by trivial algorithm, see Figure 1. 

The Algorithm 1 is not efficient for extremely symmetric matrices, such as I n : 
in that case bound step does not ever occur, because all the remaining rows 
are always equally good. Hence, Algorithm 1 must be improved, in order to 
detect some symmetries, and to avoid some unnecessary repetitions. Suppose 
that there remain / rows not included in A n , and that the column classes 
defined by S^^'" 1 ^ are such, that all column classes in the remaining rows 
are uniform (they contain either all ones or all zeros), except for at most one 
column class, which in that case has / columns, with the row and column sums 
both equal to I — 1 or 1. Then, because of the symmetry, it is enough to put 
m 

only one of the / remaining rows. After the incorporation of this 
simple heuristic, the algorithm much more efficiently deals with the matrices 
such as /„, the complement of I n , and the other highly symmetric matrices. 

Using Algorithm 1, it is possible to determine A^ for an arbitrary A e A n : it is 
enough to find ^-representatives of all (n + l) 2 matrices XiAXj, < i, j < n, 
and then to choose the smallest among them. 

One of the outputs from Algorithm 1 is the number of the pairs of row/column 
permutations, transforming A into A n . That number is used to determine the 
size of the 7r-class of A T , as it will be demonstrated below. 

Consider the problem of counting the matrices in the 7r-class of an arbitrary 
A E An- For an arbitrary B G A n let B denote the matrix obtained from B 
by sorting its rows. If A has ik groups of k equal rows, 1 < k < n, then the 
number of matrices that could be obtained from A by row permutations is 

n 

a = n\f[[ i k \ 

k=i 
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i = 0, A„ 
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= 3, £( 3 ) = {0,1,2,3}, L< 3 ) = (3) 
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Fig. 1. An example of 7r-representative determination by Algorithm 1. 

The representative of these a matrices is A . An arbitrary matrix A', obtained 
from A by a column permutation, generates in the same manner a new set of 
a matrices if and only if A' ^ Aq. If the number of different matrices A' is 
b, then the size of the 7r-class of A is ah. It is simpler to obtain b by counting 
the number p of column permutations A' of A satisfying A' = A , because 
b = n\/p. Note that p is preserved by row and column permutations of A. 

Applying Algorithm 1 to (A T ) n , p is obtained even more easily. Indeed, suppose 
that A is already a 7r-representative, i.e. A = A n . Then Algorithm 1 counts the 
row permutations A" of A, such that there exists a column permutation A'" of 
A", equal to A w . Now we find A' = ((A T ) 7r ) T and apply Algorithm 1 (again) 
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to (A') T . The matrix (A') T is a ^-representative, because ((A') T ) n = (A') T . 
Algorithm 1 gives the number of row permutations (A") T of {A') T , such that 
there exists a column permutation (A'") T of (A") T , equal to (A') T . In other 
words, we obtain the number of column permutations A" of A', such that there 
exists a row permutation A'" of A", equal to A' — which is exactly p (count 
in Algorithm 1). 

Example 6 Looking again at Example 5, we see that there are two pairs 
(P,Q) that minimize PAQ. Therefore, there are 3! 2 /2 = 18 matrices in the 
n- class of A T . 

The problem of counting the matrices in the SNF-class of an arbitrary A G A n 
is much harder. It is even harder is to enumerate the sets A n ^ = {A E A \ 
rank A = k}, < k < n: (especially m n = A n , n ) We now explicitly enumerate 
the sets A n> i, «4 n ,2, using the following characterization of matrices in A n ,2- 

Lemma 7 // the matrix A E A n> 2 contains three different nonzero columns 
a, b, c, then one of them is equal to the sum of the other two, for example 
c = a + b. Furthermore, the set of nonzero rows of the matrix [a b] equals to 
{[0 1], [1 0]}. There can not be four different nonzero columns in A. 

PROOF. Suppose A G A n ^- If two nonzero columns of A are linearly de- 
pendent, then they are obviously equal. Suppose a, b, c are the three different 
nonzero linearly dependent columns, i.e. aa + (3b + 7c = for some integers 
a, j3, 7. The coefficients a, (3, 7 must be nonzero; otherwise, if for example 
a = 0, then (3b + 7c = implies b = c. Denote by U the set of nonzero rows 
of the n x 3 matrix [a b c]. Then 

• \U\ > 1; otherwise it would be a = b = c. 

• U n {[1 0], [0 1 0], [0 1]} = 0; if, for example [1 0] G U, then a = 0. 

• therefore, U C {[1 1 1], [0 1 1], [1 1], [1 1 0]} and U ^ {[1 1 1]}. 

• [1 1 1] i U; if [1 1 1] G U, and for example [0 1 1] G U, then from a+/3+7 = 
and f3 + 7 = 0, it follows a = 0. 

• U ^ {[0 1 1], [1 1], [1 1 0]}; otherwise /3 + 7 = 0, a + 7 = 0, a + (3 = 
implies a = (3 = 7 = 0. 

Hence, there are three possibilities for U left: {[0 1 1], [1 1]}, or {[0 1 1], [1 1 0]}, 
or {[1 1], [1 1 0]}. If U = {[0 1 1], [1 1]}, then/5 + 7 = 0, a + 7 = implies 
(a, /9,7) = 7(— 1,-1,1), i.e. c = a + b; the set of nonzero rows of [a b] is 
{[0 1], [1 0]}. The two other cases are symmetrical. 

Suppose that A contains four different columns a, b, c, d. Then we must have 
for example c = a + b and the set of nonzero rows of [a b] is {[0 1], [1 0]}. 
Applying the first part of Lemma to a, b, d, we conclude that d = a + b or 
a — b + dorb — a + d. But d = a — b and d = b + a are impossible, and 
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d = a + b implies d = c. The lemma is proved. □ 

Theorem 8 a) For an arbitrary A G A n the following three statements are 
equivalent: 

(1) rank A = 1; 

(2) SNF(A) = (l.O"" 1 ); 

(3) A contains a column a ^ 0, such that all nonzero columns of A are equal 
to a. 

The number of matrices in A n> i equals 

\A n ,i\ = (2 n - l) 2 . 

b ) For an arbitrary A G A n the following three statements are equivalent: 

(1) rank A = 2; 

(2) SNF(A) = (1,1, n " 2 ); 

(3) • A contains the two nonzero columns a^b, such that all columns of A 

are in {0, a, b}, or 

• A contains the two nonzero columns a ^b, such that the set of nonzero 
rows of [a b] equals {[0 1], [1 0]} 7 and that the set of nonzero columns 
of A is {a, b, a + b}. 

The number of matrices in A n> 2 equals 

\A n>2 \ = (3" - 2 • 2 n + 1)(2 • 4" - 3 • 3 n + l)/2. 

PROOF, a) If rank A = 1 then A contains nonzero column a, such that all 
nonzero columns of A are equal to a. By subtracting one of nonzero columns 
from the others, we obtain an equivalent matrix with exactly one nonzero 
column a. By the column permutation column a is moved to the first position, 
and by the row permutation some 1 is moved to the upper left corner. By 
subtracting the first row from the other nonzero rows, we obtain that SNF of 
A is (1, 0™ -1 ). How many matrices of rank 1 there are? The number of choices 
for nonzero column a is 2 n — 1, and the number of matrices corresponding to 
the fixed a is 2™ — 1: each its column is or a, but at least one of them has to 
be equal to a. Hence, \A n ,i\ = (2 n — l) 2 . 

b) If rank A = 2 then A contains two linearly independent columns, such 
that the other columns are their linear combinations. The number of different 
nonzero columns in A is either two or it is greater than two. 

Case 1. Suppose there are exactly two different nonzero columns a, b in A. 
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The number of such matrices A is 

/9 n — 1 \ 

(3 n -2-2 n + l). 

Indeed, the number of choices for a, b equals to the above binomial coeffi- 
cient. Without loss of generality we suppose that a < b. For fixed a, 6, by 
the inclusion-exclusion principle the number of matrices A is 3™ — 2 • 2 n + 1, 

because 

• 3 n is the number of matrices with the columns from the set {0, a, b}, 

• 2 n is the number of matrices without a, and also the number of matrices 
without b, 

• 1 is the number of matrices without a and b. 

Case 2. If there are more than two different nonzero columns in A, then by 
Lemma 7 there are two different nonzero columns a, b (a < b) in A, such 
that the set of nonzero columns in A is {a, b, c = a + b}, and such that the 
row set of the matrix [a b] is {[0 1], [1 0]}. There are (3 n -2-2 n + l)/2 choices 
for columns a, b satisfying these conditions. Indeed, consider all matrices 
[a b], [b a]: 

• 3 n is the number of matrices with the row set {[0 0], [0 1], [1 0]}, 

• 2 n is the number of matrices without the row [0 1], and also the number 
of matrices without the row [1 0] , 

• 1 is the number of matrices without the rows [0 1], [1 0]). 

The number of matrices [a b] is therefore (3 n — 2 • 2 n + l)/2. The number 
4 n — 3 • 3 n + 3 • 2 n — 1 of matrices with the set of nonzero columns {a, b, c} 
(where c = a + b) is also obtained by the inclusion-exclusion principle: 

• 4™ is the number of matrices with all the columns 0, a, b, c; 

• 3 n is the number of matrices without the column a (and analogously with- 
out b, c); 

• 2 n is the number of matrices without columns a, b (and analogously with- 
out a, c; and without b, c); 

• 1 is the number of matrices without columns a, b, c. 

Therefore, the number of matrices of the rank 2, with more than two differ- 
ent nonzero columns equals 

(3" - 2 • 2 n + 1)(4™ - 3 • T + 3 • 2" - l)/2. 

The total number of matrices in A n ,i equals 

2 ( 2 ™2~ ^ + ( 4 " - 3 ' 3 ™ + 3 ' 2n - !)) (3" - 2 • 2™ + l)/2 = 
= (3" - 2 • 2™ + 1)(2 • 4 n - 3 • T + l)/2. 

In either case, in order to obtain SNF(A), the other nonzero columns are first 
transformed to by subtracting a, b or a + b from them. Next, in [a b] there is 
a row [0 1], because a < b; using that 1, the other elements of b are changed to 
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0. Finally, choosing some 1 in a, and subtracting if necessary that row from the 
others, after permuting rows/columns, we obtain the SNF. Hence, rank A = 2 
implies SNF(A) = (1, 1, n ~ 2 ). □ 

2.4 Iterative classification of (0, 1) matrices 

According to Lemma 2 we have 

Al +1 = U AeA * bord^(A). 

By changing the order of calculations, it is possible to simplify repeated de- 
termination of 7r-representatives of matrices from bord(A) by Algorithm 1. 
Matrices B in bord(A) are of the form (2). For each y the 7r-representatives 
of .B's corresponding to various inserted rows [x b] are found spending smaller 
number of steps. The point is that the rows of the 7r-representative preceding 
the row [x b] are already determined for some previous variants for that row. 

Somewhat more detailed description follows. Determine first the 7r-represen- 
tative of the matrix, corresponding to x = 0, b = 0; the inserted zero row [x b] 
is certainly the first row in the 7r-representative. The corresponding row and 
column permutations P, Q are recorded. The remaining pairs (x, b) are then 
considered in turn, lexicographically ordered. The question arises, to which 
position I might [x b] be moved during the 7r-representative determination, 
skipping the determination of first / — 1 rows of the representative. The obvi- 
ous lower bound for I is the smallest among all positions where the previous 
rows w, obtained from [x b] by changing exactly one 1 into 0, have been moved 
(except if there was an alternative to w during that step, i.e. if L^> had more 
than one member at the moment when w arrived to its destination). 

Instead of extending all A e A„, it is enough to extend the matrices from the 
set A£ of all (^-representatives in A n . By extending all A e A^ a subset of 
»4n+i is obtained; the set of ^-representatives of matrices from that subset is 
exactly A^+i- 

It is convenient to use a balanced tree to collect 7r-representatives in an or- 
dered fashion. We chose AVL tree [6] — the binary search tree satisfying the 
condition that, for every node, the difference between the heights of its left 
and right subtrees is at most 1. For n = 8, in order to save memory, a com- 
bination of AVL tree and the sorted array of matrices is used: from time to 
time the content of the tree is merged into the array. After collecting all n- 
representatives, the 7r-representatives set is reduced to the corresponding <fi- 
representatives set. To determine the set of ^-representatives, corresponding 
to a given set of 7r-representatives, the following simple algorithm is used. 
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Algorithm 2 Reduction of a given set L n of rr -representatives to the set 
of corresponding -representatives. 

{ T — auxiliary AVL tree used to collect n -representatives . } 
while L w ^ 

while there is a space in T for at least (n + l) 2 matrices 
remove the first matrix A from L n ; 

generate the set Ta of n -representatives contained in the (f>-class of A; 
insert T A into T; 
insert A^ into L^; 
remove from L n all the matrices contained in T; 

The classification of *4. 8 lasted about a month in parallel on five PC's. A huge 
number of collected 7r-representatives of order n = 8 caused serious difficulties. 
The space requirement is reduced by dividing ^-representatives into subsets, 
according to their SNF. For each extended matrix, its SNF is determined, and 
the 7r-representatives are classified into subsets with the same SNF. These sub- 
sets are then independently processed. The hardest was the SNF-class (l 7 , 0), 
with 5204144555 ^-representatives contained in a number of non disjoint sub- 
sets. These subsets were independently processed by Algorithm 2, producing 
the non disjoint sets of (^-representatives; their union consists of 71348129 0- 
representatives, approximately 1/3 of matrices in Ag. 

In order to save the space, L w and are stored in a sorted, compressed form: 
one byte for each matrix row; the group of consecutive matrices with the same 
first n — 2 rows is stored so that the common n — 2 rows are stored only once. 
As a result, the average space for a matrix of order 8 was little more than two 
bytes. 

If somebody tries to extend ^-representatives of order 8, he could expect to 
process about 300 times more ^-representatives, each giving approximately 4 
times more 7r-representatives. Therefore, the classification of matrices of order 
9 is expected to last 1000 times longer, requiring huge memory. 



2. 5 Results of classification 

We start with the simplest nontrivial case. 

Example 9 The 16 matrices of order 2 are divided into 3 <f>-classes, which 
are further subdivided into 7 n-classes: 
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In Table A.l all the 36 ir- representatives of order 3 are shown. The 5 SNF- 
classes are in separate blocks, divided into compartments with ^-classes. The 
first matrix in each 0-class is the smallest 7r-representative, i.e. the (^-rep- 
resentative. For each 7r-and SNF-class, their size is given. The matrices are 
represented by hexadecimal vectors, each component representing a row of a 
matrix. For example, the last vector (3, 5, 6) in Table A.l represents the matrix 

1 1 
10 1 • 

1 1 

The matrix (1,2,5) is a ^-representative of the matrix from Example 5. 

In Table A. 2 all the 39 ^-representatives of order 4 are shown, together with 
the sizes of their 0-classes. 

In Table 2 p n , \A^\, s n , a n , \T> n \, and the set V n are given for 1 < n < 8, 
where s n = \S n \ and p n = (2 n /(n + 1)! 2 )/|.4*|. In the last row of Table 2 s 9 , 
\T>g\, a$, Vg are given; the explanation how they are obtained will be given in 
section 3. 

Denote by F{n) the following statement: 



A E A n , satisfying SNF(A) = d = (di, d 2 , . . . , d n ) exists if and only if(6) 
there exists A' e A n +i, satisfying SNF(A') = d' — (d± : d 2 , ■ ■ ■ , d n , 0). 

Obviously, the first condition implies the second one. The implication in the 
opposite direction is not obvious at all; it would follow from the following 
stronger statement: 



H(n) : Let A' e A n+1 , rank A' = n, and SNF(A') = d' = (d u d 2 , . . . , d n , 0). 
Then A' has at least one minor A E A n with SNF(A) = d = (d 1: d 2 , ■ ■ ■ , d n ). 
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Table 2 

The numbers of equivalence classes in A, 



n 


Pn 


1 A 4* 1 
\*A-n \ 




a n 


\D n \ 


v n 


1 


0.250 


2 


2 


2 


2 


{0,1} 


2 


0.148 


3 


3 


2 


2 


{0,1} 


3 


0.074 


12 


5 


3 


3 


{0-2} 


4 


0.117 


39 


8 


4 


4 


{0-3} 


5 


0.167 


388 


14 


6 


6 


{0-5} 


6 


0.334 


8102 


26 


10 


10 


{0-9} 


7 


0.528 


656103 


56 


19 


22 


{0 - 18,20,24,32} 


8 


0.701 


199727714 


129 


41 


46 


{0-40,42,44,45,48,56} 


9 






333 


103 


114 


{0 - 102,104,105,108,110, 
112, 116, 117, 120, 125, 128, 144} 



But the following matrix F e A\q is a counterexample to if (10): 



11 


10 10 1 


110 


10 10 1 


110 


110 10 


10 1 


10 110 


11 


111 


110 


1110 


10 10 


1110 


10 1 


1110 


110 


110 1 


10 1 


10 11 



The matrix F consists of blocks A, B, C, D, having 2, 3, 2, 3 ones in each row 
respectively, and also having 2, 2, 3, 3 ones in each column, respectively; F is 
singular, because the sums of rows of [A B] and [C D] are equal. It can be 
verified that rankF = 9, SNF(F) = (1 9 ,0), but all minors of F have SNF 
different from (l 9 ). 

In Table A. 3 the SNF-representatives of matrices in A n , n < 8, are listed, ac- 
companied with the size measures of corresponding SNF-classes (the number of 
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Table 3 



The number of matrices of the rank k in A n , n < 8. 



n 

k 


1 2 


3 4 5 6 


7 


8 





1 1 


111 1 


1 


1 


1 


1 9 


49 225 961 3969 


16129 


65025 


2 


6 


288 6750 118800 1807806 


25316928 


336954750 


3 




174 36000 3159750 190071000 


9271660734 


397046059200 


4 




22560 17760600 5295204600 


1001080231200 


144998212423680 


5 




12514320 34395777360 


32307576315840 


17952208799918400 


6 




28836612000 259286329895040 


720988662376725120 


7 






270345669985440 


7547198043595392000 


8 








10160459763342013440 



matrices, the number of 7r-representatives and the number of ^-representatives 
in each SNF-class). The sizes of 7r-classes are determined using Algorithm 1. 
The classes are ordered lexicographically by the SNF (with zeros moved to 
the end of SNF). 

One can verify this classification starting from the sorted list of all (^-repre- 
sentatives. For each of them one has to check if it is indeed a ^-representative. 
The next step is to sum the numbers of 7r-representatives in all 0-classes, and 
to compare the sum with the corresponding entry in Table 1. One could also 
check that the sum of sizes of SNF-classes in A n equals 2 n for each n < 8, see 
Table A. 3. The sorted lists of ^-representatives for n < 8 can be downloaded 



from http: //www.matf.bg.ac.yu~ezivkovm/01matrices.htm . 



We now review some interesting facts, which are seen from Table A. 3. 

Let T(n,k) = \A n ^\- In Table 3 the numbers T(n,k), < k < n < 8, 
are shown (of course, they are easily obtained from Table A. 3). The part of 
Table 3 corresponding to n < 7 is the same as in [10]; it is also an entry 
in [4, Sequence A064230]. Another interesting entry in [4, Sequence A055165] 
is the sequence m n , where m n is the number of regular (0, 1) matrices of order 
n — the diagonal of Table 3. The seemingly new member of that sequence 
is m 8 = 10160459763342013440. If we suppose that all matrices in A n are 
equiprobable, then the rank probability distribution is shown in Table 4 for 
n < 8. Looking at Table 4, one could erroneously conclude that large fraction 
of matrices in A n is singular. In fact, the fraction of singular matrices in A n 
tends to for n large [11]. 
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Table 4 



The probability that a random matrix in A n has the rank k, < k < n < 8. 



n 


1 


2 


3 


4 


5 


6 


7 


8 


k 





















0.5 


0.0625 


0.00195 


0.00002 


0.00000 


0.00000 


0.00000 


0.00000 


1 


0.5 


0.5625 


0.09570 


0.00343 


0.00003 


0.00000 


0.00000 


0.00000 


2 




0.3750 


0.56250 


0.10300 


0.00354 


0.00003 


0.00000 


0.00000 


3 






0.33984 


0.54932 


0.09417 


0.00277 


0.00002 


0.00000 


4 








0.34424 


0.52931 


0.07706 


0.00178 


0.00001 


5 










0.37296 


0.50052 


0.05739 


0.00097 


6 












0.41963 


0.46059 


0.03908 


7 

8 














0.48023 


0.40913 
0.55080 



Table 5 

The possible numbers of 7r-orbits inside ^-orbits of A n - 



n 


The set of </>-orbit sizes 


1 


{1} 


2 


{1,2,4} 


3 


{1,2,4,5,9} 


4 


{1-5,7,9-11,13,16,17} 


5 


{1 - 18,20,21,25,26,30,36} 


6 


{1, 2, 4 - 27, 29 - 32, 35 - 37, 42, 49} 


7 


{1 - 38, 40, 42 - 44, 48 - 50, 56, 64} 


8 


{1 - 46, 48 - 51, 53, 54, 56 - 58, 63 - 65, 72, 81} 



It turns out that F(n) (6) is true for n < 7, i.e. the set of SNF's of rank k 
is the same for all n, k < n < 8. For example, the SNF-representative of the 
SNF-class (1, 1, 2, n " 3 ) is the matrix (0 n " 3 , 3, 5, 6) for 3 < n < 8. 

The smallest n for which there are two matrices in A n with the same deter- 
minant, but with different SNF's is 5: SNF(3,C,15,16,19) = (1,1,1,4) and 
SNF(3,5,9,11,1E) = (1,1,2,2). 

In Table 5 the possible numbers of 7r-orbits inside 0-orbits are shown for 
1 < n < 8. These numbers are between 1 and (n + l) 2 ; as it is seen, the value 
(n + l) 2 is attained only if n > 5. 
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Table 6 

The maximal ADV's of matrices from A n +i, obtained by extending matrices equiv- 
alent to I n . 



n 


| detA\ 


A 


3 


3 


3 5 9 E 


4 


5 


3 5 E 16 19 


5 


9 


3 D 15 IA 26 39 


6 


18 


7 19 2A 34 AC 53 65 


7 


40 


7 19 56 65 9C S3 


8 


105 


7 39 5A AC D5 E3 136 UD 19 B 



If A G A n , A ~ I n and 5 G bord(A), then SNF(£?) contains at least n ones, 
see (5). The question arises, what are the possible values of the last element 
of SNF(_B), i.e. which values can take | detS|? The largest possible values of 
| det.B| under these assumptions, along with the examples of matrices B for 
which these values are attained, are given in Table 6. In fact, the matrices 
from Table 6 maximize | det Bj det A\ for all regular A G A, n < 8. 

More generally, it is interesting to describe the relationship of SNF(A) to 
SNF(A') if A' G bord(A). During iterative classification, the sets 

{SNF(5) | B G bord(A), A G A n , SNF(A) = s} 

are recorded for all SNF-classes s G 5„. The results are represented by the 
incidence matrix M n of dimensions |«S n | x |«S n +i|, with entries 

1, if there exist A G A n and B G bord(A), with SNF's s and s' 

0, otherwise 

(7) 

Let G(n), denote the following statement: 



G(n) : There exist matrices A G A n , A' G bord(A), such that (8) 
SNF(A) = (d 1 ,d 2 ,..., d n ), SNF(A') = (d[, d' 2 , . . . , d' n , d' n+1 ) if and only if 
there exist matrices B G »4 n+ i, B' G bord(5), such that 
SNF(5) = (d ± , d 2 , . . . , ck, 0), SNF(B') = (d[, d' 2 , ...,d' n , d' n+1 , 0). 

By exhaustive search it is verified that G(n) is true for n < 6, enabling to put 
all the transposed incidence matrices M n , n < 7 together into single Table A. 4. 
The l's are represented by • ; the 0's are represented by * if they are the 
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consequence of the following Lemma (describing constraints for SNF(A') if 
A' E bord(A)); otherwise, they are represented by o . 



Lemma 10 For an arbitrary A E A n , let A' E bord(A), and let SNF(A) = 
(di, d 2 , . . . , d n ), SNF(A') = «, d' 2 , ...,d' n , d' n+1 ). Then 

(1) rank A < rank A' < rank A + 2; 

(2) d! x d' 2 . . . d\ divides d\d 2 . . .d; L for all i, 1 < i < rank A; 

(3) IVl=idi divides det A'. 



PROOF. 

(1) The first inequality follows from the fact that the rank of a sub matrix 
is a lower bound on the rank of a matrix. The second inequality follows 
from the observation that A' is an at most rank 2 perturbation of A. 

(2) This is a direct consequence of the fact that d\d! 2 . . . d\ is the largest 
common divisor of all minors of A' of order i, see for example [7]. 

(3) Let P, Q be the matrices such that SNF(A) = PAQ = D = d 2 , . . . , d n ), 
| detP| = | det (5| = I. Let 

^ y 

A' = 

x b 

The case det A' = is trivial; suppose det A' ^ 0. If xQ = [ai a 2 ... a n ], 
Py — [ci c 2 ... c n ] T , then from the identity 



P 




Ay 




Q 




DPy 


1 




x b 




1 




xQ b 



it follows (another way to express determinants of matrices obtained by 
extension, see (3)) 

n n 

det A' = b]ld i -Y i a i c i J] d r ( 9 ) 

i=l 1 = 1 l<j<n, 

Since rank A' — n + 1, then we have rank A > n — 1. If rank A = n — 1, 
then d n — 0, implying 

n-l 

det A' = -a„c„ ]J djj 
i=i 

otherwise 

(n \ n—1 

bd n - CLiddn/ di j JJ dj. 
i=i / i=i 
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In both cases niLi 1 di divides det A. □ 

Suppose A G An, A' G bord(A). From Table A. 4, we see the following inter- 
esting facts: 

• The first o in some M n corresponds to s = (1, 0), s' — (1, 1, 2). It is equiva- 
lent to following statement: if A G A2,i then | det A'\ < 2. 

• if A G A 3 , 2 , then | det A'\ < 3. 

• if A G At', SNF(A) = (1,1, 1, 0), then | det A'\ < 5. 

• if A G Aa, SNF(A) = (1, 1, 2, 0), then | det A'\ < 4. 

• if A G A b , SNF(A) = (1,1, 1, 1, 0), then | det A'\ ^ 7. 

• if A G A 5 , SNF(A) = (1,1, 1, 2, 0), then | det A'\ ^ 6. 

• if A G A, SNF(A) = (1,1, 1, 3, 0), then | det A'\ £ {6, 9}. 

• if A G A, SNF(A) = (1,1, 1, 2, 2), then SNF(A') ^ (1, 1,1,1, 4, 0). 

• if SNF(A) = (l"" 1 ,^) and SNF(A') = (l"" 1 ,^^) then d' n divides d n for 
all n < 7. 

• if SNF(A) = s = (l" -1 ,^) G «S n and SNF(A') = s' = (l n ,d' n+1 ) G S n+1 
then 

• if n < 6, then m SiS / = 1. 

• if n — 7, then m S)S / = 1 if and only if 

(d n , C) ^ {( 17 - 34 )' ( 7 - 39 ). ( 13 - 39 ). (!» 42 )' 

(4, 42), (6, 42), (7, 42), (13, 42), (14, 42)}. 

• if n — 8, then there are more exceptions to m S;S ' = 1, but there is one 
exotic group of them: if d n = 19 then d' n+1 must be divisible by 19; 19 is 
the only integer satisfying such a condition. 



3 Determinant and SNF sets of (0, 1) matrices of order 9 

Determination of {| det(A')| | A' G bord(A)} is a simple operation, see the ex- 
planation following (3). It was effectively performed for all 199727714 matrices 
in Ap, merging these sets V 9 is obtained, see Table 2. 

The similar idea — determine ADV's, and only if necessary, determine SNF's 
of the results of extension — is used to obtain <S 9 . Suppose we know in advance 
the number f d of different SNF's in V 9 corresponding to a given ADV d > 0. 
During the extension of matrices from A n , the SNF's of extended matrices 
with the ADV d are determined only if the number of SNF's with ADV d 
is still less than fy. If we know only upper bound on fd, then the heuristic 
does not work — we have to determine SNF's of all matrices with the ADV 
d. Therefore, it is useful to determine f d for at least some d > 0. 
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Table 7 

The number of partitions of r into at most n positive integers. 



r 


12 3 4 5 6 7 


n 







10 


1 


111111 1 1 


2 


112 2 3 3 4 4 


3 


112345 7 8 


4 


1 1 2 3 5 6 9 11 


5 


1 1 2 3 5 7 10 13 


6 


1 1 2 3 5 7 11 14 


7 


1 1 2 3 5 7 11 15 



Denote by p n {r) the number of partitions of r into at most n positive integers. 
In order to determine the upper bound for fy, suppose first that d is a prime 
power, d = p r . If A e ^4 n and | det A\ = d, then SNF(A) is of the form 

n 

(p Xl ,p X2 ,...,p Xn ), 0<x 1 <x 2 <---<x n , J2 x i = r - 

1=1 

The number of different exponent vectors (x±, x 2 , ■ ■ ■ , x n ) is equal to p n (r). 
The values p n (m) are computed using the recurrence (see for example [12]) 
p n (0) = 1, n > 0, p (r) = for r > 1, and p„(r) = p n -i(r) + p n (r - n), see 
Table 7. 

Example 11 // d = 8 = 2 3 and n = 6 we nave pe(3) = 3; SNF(A) is one of 
(1 5 ,8) 7 (1 4 ,2,4) and (l 3 , 2, 2, 2). We see /rom Table A. 3 that all these SNF do 
exist, i.e. for each of them there exists some (0,1) matrix. Another example 
d = 3 2 , n = 6, shows that pe(2) = 2 is only an upper bound: the SNF-class 
(l 4 , 3, 3) is empty. 

More generally, if d = Il^liPfS where are different primes, then the upper 
bound on the number of different SNF's with the ADV d is Y\ZLiPn{oti). 

Example 12 If n = 8 and d = 36, then there are p$(2)p 8 (2) = 4 such SNF's: 
(1 6 ,2, 18), (l 7 , 36), (1 6 ,3,12), (l 6 ,6,6);a// these SNF's are found in Table A. 3. 

In order to obtain a tighter upper bound for the number of different SNF's, 
we have to include somewhat more information. If we further suppose that 
A' G bord(A) and SNF(A) = s = (di, d 2 , ■ ■ ■ , d n ), then by Lemma 10 for some 
s' = (d[, d' 2 , . . . , d' n+1 ) the equality SNF(A') = s is impossible. For example, if 
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s contains k ones, then s' contains at least k ones and if rank A > k + 1 then 
d' k+1 divides d k+1 . 

Using these facts, the regular part of Sg was determined, see Table A. 5. The 
^-representatives from the chosen SNF-class of As were extended computing 
determinants, and, if necessary, determining SNF's. The upper bounds for the 
number of different SNF's, obtained by Lemma 10, are rough for larger ADV 
values, but the consequences are not dangerous, because of the small number 
of extended matrices with the large ADV: it is not hard to compute the SNF's 
of all of them. 

To complete <Sg, it is necessary to determine the singular part of Sg. If we 
would know that F(8) is true, then the set of singular SNF's of order 8 would 
be equal to S$ (with each SNF extended by one zero, of course). Not knowing 
a simple proof of F(8), we proceed with a shortened exhaustive proof. 

The idea is to narrow the set of SNF-classes in <S 8 , the extension of which 
can lead to a new singular SNF of order 9. If SNF(A) = d 2 , ■ ■ • , d 8: 0) 
for some A e Ag, then (because we know the set of SNF's of lower or- 
ders) by Lemma 10 we can narrow the set SNF-classes, containing A. We 
obtain that the only new possible SNF's are the following SNF's of the rank 8: 
(l 7 , m, 0), m = 44, 45, 48, 56 and (l 6 , 2, 28, 0); and the following SNF's of the 
rank 7: (l 6 , 20, 0,0), (l 6 , 24, 0,0), (l 6 , 32, 0,0), (l 5 , 2, 12, 0, 0), (l 5 , 2, 16, 0, 0), 
(l 5 , 4, 8, 0,0), (l 4 , 2, 2, 8, 0,0), (l 4 , 2, 4, 4, 0, 0). 

The extension of which matrices gives the matrices with such SNF's? For ex- 
ample, we know that the SNF (l 7 , 44, 0) can be obtained only by the extension 
of a matrix in which 44 divides all minors of order 8; therefore 44 also divides 
a nonsingular minor of order 8; hence the SNF of that minor could be only 
(l 6 , 2, 22). Considering analogously the rest of listed SNF's of order 8, we ob- 
tain that matrices from Ag,s, with the SNF equal to some from the list above, 
can be obtained only by the extension of matrices from As with the SNF 
(l 6 , 2, 22), (l 6 , 2, 24), (l 6 , 3, 15), (l 5 , 2, 2, 12), or (l 5 , 2, 2, 14). 

Analogously, we obtain that matrices from Agj with one of the listed SNF's, 
can be obtained only by double extension of matrices from A7 with the SNF 
(1 5 ,2, 10), (1 4 ,2,2,6), or (l 3 , 2, 2, 2, 4). After the complete search through all 
matrices that can be obtained by the extensions listed, it is found that there 
are no new singular SNF's of order 9 i.e. that F(8) is also true. That completes 
the determination of Sg. 

In Table A. 6 the part of the incidence matrix M§ is shown, corresponding to 
regular matrices in <S 9 . The table was obtained by extending ^-representatives 
from As j and ^4 8 ,8 ; the singular extended matrices were ignored. 
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4 The lower bounds for the first missing determinant, a, 



Denote by f n the nth Fibonacci number (f\ — f 2 — 1 and f n = / n _i + / n _ 2 
for n > 3). Paseman [13] shows that a n > 2f n _ ± . We give the sketch of his 
proof, and then we give the sharper lower bounds for a n , n < 19. 

Consider the so called Fibonacci matrices F n G A n with the element 
equal to 1 if and only if j — % = —1,0,2,4, . . .; det F n = f n . The cofactors 
corresponding to the first row of F n are / n _i, /„_ 2 , -f n -3, ~fn-4---,-fi- 
Consider the matrix U G bord(F n ), 



Fn y 
x b 

where x = \x x x 2 ■■■ x„_i], y = [yi y 2 • • • y«-i] T - Let yi = 1, y 2 = yz = • • • = 
y n = and rri = x 2 = 0. Then from (3) we have 

n-2 

det U = ^ x n+ i_i/i + &/„. 
i=i 

Therefore, each integer from [0, 2f n — 1] is determinant of some U G bord(F n ), 
and a n > 2/ n _i. 

In order to prove that a n > m, one can give a list of matrices from A n -i, such 
that determinants of their extensions cover [l,m — 1]. The proof verification 
then includes the procedure of finding determinants of all extensions of a given 
matrix. Still, such a list is essentially more compact than the list of matrices 
from An, with determinants covering [l,m — 1]. 

Denote by the minimal integer not in U{|det5| | B G bord(A)}, the 
"extension spectrum" of A G A n . In this context, the matrices A with high 
cia are of special interest. If > 1 and SNF(A) = (di, d 2 , . . . , d n ), then 
c?i = d 2 = ■ ■ ■ = d n _i = 1, because determinants of all extensions of A are 
divisible by d n _i, see (9). 

In order to find lower bounds for some a n , one can start from a well chosen 
set B n -i C A n -i, and then to find ADV's of all extended matrices. If m is the 
smallest number not equal to some of these ADV's, then a n > m. Afterwards, 
some subset of extended matrices with different SNF's is taken to be the set 
B n , and the next iteration can be started. 

The starting set Bg was constructed in the following way. From each SNF-class 
in ^8 a number of matrices is taken, with different numbers of 7r represen- 
tatives in their 0-classes. Extending these matrices, a set of matrices with 
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different SNF's is obtained, but without any matrix with the SNF (1 8 ,97). 
By adding one such matrix, the set Bg is completed. The sets B±o, Bu an d 
£>i2 are generated iteratively, as explained above. At the end, the ADV's of 
all matrices obtained by extending the matrices in B 12 are determined. The 
resulting lower bounds are a w > 259, an > 739, ai 2 > 2107, a 13 > 6157. 

For n > 13 we used an alternative heuristic, described by Algorithm 3. 

Algorithm 3 Heuristic to find lower bound for a n+1 . 

Input : C n C A n , list of matrices to be extended. 
Output : lower bound for a n+ i, and list C n+ i C A n+ i 
of "promising" matrices for the following iteration. 
{ Initialization: } 

firstO <— 1; { the first integer not "covered" by ADV's } 
dmax <— 1; { the largest AD V found until now } 
£ n +i <— 0/ { output list } 
for all A G C n do 



{ Consider the extensions A' 



} 



A y 
x b 

Compute det A and B = [B^] = adj A; { transposed cof actor matrix of A} 
for all y E {0, 1}" do 

{ the next linear combination of rows of B } 

determine the coefficients of the linear combination 

-bdetA + TZ =lXi (v; ;//,/;,,) 

and the sums s + , s~ of its positive and negative members; 
if max{— s~,s + } > firstO then {"poor" linear combinations are skipped} 
for all (x,b) E {0, l} n+1 do 

compute det A' ; { by one addition only, using Gray code } 
if | det A'\ = firstO then 

update firstO; 
if | det A' | > dmax then 

dmax <— | det A' | ; 
if | det A'\ > 0.9 dmax then 
append A' to C n+ \; 



Elimination of "poor" linear combinations is a powerful heuristics if the ma- 
trices with the high extension spectra are placed in the beginning of C n . The 
major part of linear combinations is skipped after only a few first matrices in 
C n , reducing the extension complexity roughly to 0{n2 n ) (instead of 0(4 n )). 
In Table A. 7 for 10 < n < 19 we give 

• lower bound for a„, 
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• IAj-iI, the number of extended matrices, 

• a matrix A n _i with the highest extension spectrum found in |^4 n _ 1 |, 

• extension spectrum and determinant of A n _i. 

Complete lists of matrices, whose extension determinants prove these lower 



bounds, can be fount at http:/ /www.matf.bg.ac.yu~ezivkovm/01matrices.htm 



5 Counting (0, 1) matrices with the maximum determinant 



Using the classification of A n , it is not hard to compute the number c n [4, 
Sequences A051752] of matrices in A n with the maximal determinant d n (i.e. 
1/2 of the number of matrices with the ADV d n ) for n < 9. 

The first 8 members of the sequence c n are found in Table A. 3; the number 
c 8 = 195955200 is new. 

In order to determine eg, from Table A. 4 we see that the matrix from Ag with 
the ADV 144 could be obtained only by extending matrices from As with the 
SNF (l 5 , 2, 2, 6) or (l 5 , 2, 2, 12). After the extension of these two SNF-classes, 
it turned out that there is a unique 0-class with the ADV 144 — the class with 
the representative (F,33,C3,FC,155,15A,166,196,1A9). Half of the number of 
matrices in that 0-class is eg = 13716864000. It is interesting that for all n < 9 
there is a unique 0-class with the maximal ADV. 
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Table A.l 

7r-representatives of 



SNF 



size 
1 





1 





SNF 




size 


1 








49 








1 


9 








7 


3 


1 


1 


1 


3 


7 


7 


7 


1 








3 


9 


3 


3 


3 


3 





1 


1 


9 





7 


7 


3 





3 


3 


9 



, 1) matrices of order 





SNF 




size 


1 


1 





288 





1 


2 


18 





1 


7 


18 


1 


1 


3 


18 


1 


1 


6 


9 


3 


7 


7 


9 





1 


3 


36 





1 


6 


18 





3 


7 


18 


1 


1 


2 


18 


1 


3 


3 


18 


1 


1 


7 


9 


3 


3 


7 


9 


1 


6 


6 


9 


1 


7 


7 


9 





3 


5 


18 


3 


3 


5 


18 


1 


2 


3 


18 


1 


6 


7 


18 



3. 





SNF 




size 


1 


1 


1 


168 


1 


2 


4 


6 


1 


2 


7 


18 


1 


3 


5 


18 


1 


3 


6 


36 


3 


5 


7 


18 


1 


2 


5 


36 


1 


3 


7 


36 



SNF 
1 1 2 


size 
6 


3 5 6 


6 



Table A.2 

(^-representatives of (0, 



1) matrices of order 4. 





SNF 




size 














1 














1 




SNF 




size 


1 











225 











1 


25 











3 


50 








1 


1 


50 








3 


3 


100 




SNF 




size 


1 


1 








6750 








1 


2 


200 








1 


3 


400 








1 


6 


600 








3 


5 


600 








3 


7 


300 





1 


1 


2 


600 





1 


1 


6 


450 





1 


2 


3 


600 





1 


3 


3 


300 





1 


6 


6 


900 





1 


6 


7 


900 





3 


3 


5 


900 





SNF 




S1ZG 


1 


1 


1 


o 


35400 





1 


2 


4 


600 


A 




1 


2 


5 


3600 





1 


2 


1-7 

7 


1800 





1 


3 


r 




1800 


n 


1 


3 


a 

yj 


36DD 
ouuu 





1 


3 


7 


3600 





1 


6 


A 


3600 





3 


5 


7 


1800 





3 


5 


9 


1200 





3 


5 


A 


3600 


1 


2 


3 


4 


3600 


1 


2 


4 


7 


1200 


1 


2 


5 


6 


3600 


1 


6 


7 


A 


1800 



SNF 
112 


size 
600 


3 5 6 


600 





SNF 




size 


1 


1 


1 


1 


20040 


1 


2 


4 


8 


600 


1 


2 


4 


9 


7200 


1 


2 


5 


A 


1440 


1 


2 


5 


B 


7200 


1 


2 


7 


B 


3600 



SNF 
1112 


size 
2400 


1 6 A C 


2400 



SNF 
1113 


size 
120 


3 5 9 E 


120 
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Table A.3 

The representatives and the sizes of SNF-classes in A n , n < 8. 





A 


l 


The number of 


The SNF-class 




det 


SNF 


matrices 


7r-classes 


(^-classes 


representative 











1 


1 


1 





1 


1 


1 


1 


1 


1 


1 


Total: 


2 


2 


2 








A 


2 


The number 


of 


The SNF-class 




det 


SNF 


matrices 


7r-classes 


(^-classes 


representative 











1 


1 


1 





1 





1 


9 


4 


1 


1 


2 


1 


1 1 


6 


2 


1 


1 2 


Total: 


16 


7 


3 







A 


3 


The number of 


The SNF-class 




det 


SNF 


matrices 


7r-classes 


^-classes 


representative 








00 


1 


1 


1 


00 


1 





00 1 


49 


9 


4 


00 1 


2 





01 1 


288 


18 


4 


01 2 


3 


1 


11 1 


168 


7 


2 


12 4 


4 


2 


11 2 


6 


1 


1 


35 6 


Total: 


512 


36 


12 





Aa 


The number 


of 


The SNF-class 




det 


SNF 


matrices 


7r-classes 


</>-classes 


representative 








0000 


1 


1 


1 





1 





0001 


225 


16 


4 


1 


2 





0011 


6750 


84 


12 


1 2 


3 





0111 


35400 


150 


14 


12 4 


4 





0112 


600 


5 


1 


3 5 6 


5 


1 


1111 


20040 


49 


5 


12 4 8 


6 


2 


1112 


2400 


10 


1 


1 6 A C 


7 


3 


1113 


120 


2 


1 


3 5 9 E 


Total: 


65536 


317 


39 
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Table A.3 
Continued 



A 


The number of 


The SNF-class 




uet 




matrices 


7r-classes 


0-classes 


representative 








00000 


1 


1 


1 














1 





00001 


961 


25 


9 











1 


2 





0001 1 


118800 


260 


37 








1 


2 


3 





00111 


3134400 


1346 


113 








1 2 


4 


4 





0112 


25350 


25 


5 








3 5 


6 


5 





01111 


16853400 


2589 


141 





1 


2 4 


8 


6 





01112 


880200 


210 


17 





1 


6 A 


C 


7 





01113 


27000 


15 


2 





3 


5 9 


E 


8 


1 


11111 


9702720 


831 


39 


1 


2 


4 8 


10 


9 


2 


11112 


2427840 


254 


15 


1 


2 


C 14 


18 


10 


3 


11113 


289440 


51 


5 


1 


6 


A 12 


1C 


11 


4 


11114 


65520 


12 


2 


3 


5 


9 11 


IE 


12 


5 


11115 


7200 


3 


1 


3 


5 


E 16 


19 


13 


4 


11122 


21600 


2 


1 


3 


C 15 16 


19 


Total: 


33554432 


5624 


388 





A 6 


The number ot 


The SNF-class 




det 


SMF 


matrices 


7r-classes 


<p-classes 


representative 








000000 


1 


1 


1 














1 





000001 


3969 


36 


9 











1 


2 





000011 


1807806 


660 


76 











1 2 


3 





000111 


189336000 


7586 


472 











1 2 4 


4 





000112 


735000 


86 


10 











3 5 6 


5 





001111 


5168108400 


47605 


1913 








1 


2 4 8 


6 





001112 


124744200 


2120 


115 








1 


6 A C 


7 





001113 


2352000 


91 


9 








3 


5 9 E 


8 





011111 


30991962960 


112080 


3262 





1 


2 


4 8 10 


9 





011112 


3122915040 


14986 


511 





1 


2 


C 14 18 


10 





011113 


226603440 


1618 


75 





1 


6 


A 12 1C 


11 





011114 


38419920 


307 


16 





3 


5 


9 11 IE 


12 





011115 


3175200 


46 


3 





3 


5 


E 16 19 


13 





011122 


12700800 


78 


4 





3 


C 


15 16 19 


14 


1 


linn 


18480102480 


39637 


952 


1 


2 


4 


8 10 20 


15 


2 


111112 


7737327360 


17642 


442 


1 


2 


4 


18 28 30 


16 


3 


111113 


1537446960 


4079 


128 


1 


2 


C 


14 24 38 


17 


4 


111114 


628548480 


1685 


52 


1 


6 


A 12 22 3C 


18 


5 


111115 


127224720 


429 


18 


1 


6 


A 1C 2C 32 


19 


6 


111116 


93139200 


263 


9 


3 


5 


9 


16 2E 31 


20 


7 


111117 


12877200 


54 


3 


3 


5 


9 


IE 2E 31 


21 


8 


111118 


6703200 


27 


2 


3 


5 


E 


19 29 36 


22 


9 


111119 


1058400 


7 


1 


3 D 15 1A 26 39 


23 


4 


111122 


208857600 


473 


17 


1 


6 18 2A2C 32 


24 


8 


111124 


3175200 


12 


1 


3 


C 15 1A 26 29 


25 


8 


111222 


151200 


2 


1 


7 19 IE 2A 2D 33 


Total: 


68719476736 


251610 


8102 
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Table A.3 
Continued 



A 7 


The number of 






The SNF-class 




det 


SNF 


m atri ppc 


7T-p1 a ccpc 






representative 








n n n n n n n 

0000000 


l 


1 


1 


U 














1 





r\ r\ r\ r\ r\ r\ ~i 

0000001 


16129 


49 


16 














1 


2 





0000011 


25316928 


1428 


170 














12 


3 





00001 1 1 


9254300328 


31994 


1908 














12 4 


4 





00001 12 


17360406 


246 


34 














3 5 6 


5 





0001111 


989588124000 


501563 


17596 











1 


2 4 8 


6 





0001112 


11359807200 


13645 


694 











1 


6 A C 


7 





0001113 


132300000 


400 


30 











3 


5 9 E 


8 





0011111 


30826279895040 


4358421 


105808 








1 


2 


4 8 10 


Q 


n 
u 


Ofll 1 1 1 9 




^1 fiQflzL 
oioyu^t 


9295 








1 


2 


C 14 18 


10 





0011113 


64153434240 


22902 


ooo 








1 


6 


A 12 1C 


11 





0011114 


8175222720 


3714 


168 








3 


5 


9 11 IE 


12 





0011115 


509443200 


413 


23 








3 


5 


E 16 19 


13 





0011122 


2694686400 


1032 


58 








3 


C 


15 16 19 


14 





0111111 


219225571810560 


13834240 


261882 





1 


2 


4 


8 10 20 


15 





0111112 


34159997168640 


2624469 


53874 





1 


2 


4 18 28 30 


16 





0111113 


4018162256640 


376699 


8633 





1 


2 


C 


14 24 38 


17 





0111114 


1176364465920 


123510 


3024 





1 


6 


A 12 22 3C 


18 





0111115 


182858215680 


23489 


633 





1 


6 


A 1C 2C 32 


19 





0111116 


110954188800 


13823 


361 





3 


5 


9 


16 2E 31 


20 





0111117 


12940704000 


2133 


64 





3 


5 


9 IE 2E 31 


21 





0111118 


5966553600 


1006 


33 





3 


5 


E 19 29 36 


22 





0111119 


829785600 


189 


7 


3 D 15 1A 26 39 


23 





0111122 


389700057600 


37489 


927 





1 


6 18 2A2C 32 


24 





0111124 


2857680000 


415 


19 





3 


C 15 1A 26 29 


25 





0111222 


127008000 


29 


4 





7 19 IE 2A 2D 33 



34 



Table A.3 
Continued 



A 7 


The number of 




The SNF-class 
representative 




det 


SNF 


matrices 


7r-classes 


^-classes 


26 


1 


linn i 


135491563468800 


5593528 


91764 


1 


2 4 8 10 20 40 


27 


2 


linn 2 


83220427382400 


3493129 


58179 


1 


2 4 8 30 50 60 


28 


3 


mm 3 


23436399974400 


1020752 


17707 


1 


2 4 18 28 48 70 


29 


4 


mm 4 


13285672243200 


581948 


10189 


1 


2 C 14 24 44 78 


30 


5 


mm 5 


3754520017920 


172714 


3169 


1 


2 C 14 38 58 64 


31 


6 


mm 6 


4201407745920 


185688 


3320 


1 


6 A 12 2C5C 62 


32 


7 


mm 7 


813250851840 


39068 


749 


1 


6 A 12 3C 5C 62 


33 


8 


mm 8 


693389168640 


32490 


645 


1 


6 A1C 32 52 6C 


34 


9 


mm 9 


257766405120 


12609 


253 


1 


61A2A 34 4C 72 


35 


10 


111 111 10 


215881142400 


10094 


199 


3 


5 9 IE 2E 4E 71 


36 


11 


liiimi 


49798425600 


2598 


55 


3 


5 9 IE 31 51 6E 


37 


12 


111 111 12 


67511808000 


3263 


71 


3 


5 E 16 39 59 66 


38 


13 


1 1 1 111 13 


12283084800 


686 


17 


3 


5 E 19 36 56 69 


39 


14 


111 111 14 


12260505600 


615 


12 


3 


5 19 29 36 4E 71 


40 


15 


11111115 


4064256000 


215 


6 


3 


D 15 26 38 5E 61 


41 


16 


11111116 


2235340800 


143 


6 


3 


C 15 36 39 5A 65 


42 


17 


11111117 


406425600 


27 


2 


3 


D 16 2E 39 5A 65 


43 


18 


111 111 18 


541900800 


24 


1 


7 


19 2A 34 4C 53 65 


44 


4 


111112 2 


4413330432000 


184475 


3220 


1 


2 C 30 54 58 64 


45 


8 


111112 4 


343226419200 


15119 


317 


1 


6 18 2A 34 4C 52 


46 


12 


111112 6 


21946982400 


997 


28 


3 


5 19 29 36 4E 51 


47 


16 


111112 8 


1219276800 


102 


6 


3 


C 31 55 5A 66 69 


48 


20 


111 1 1 2 10 


135475200 


10 


1 


7 


19 2A 34 4C 52 63 


49 


9 


111113 3 


29940019200 


1358 


37 


3 


5 18 28 49 4E 71 


50 


18 


111113 6 


139708800 


10 


2 


3 ID 2D 36 3A 4E 71 


51 


16 


111114 4 


254016000 


19 


3 


3 


C 35 3A 55 66 69 


52 


8 


111122 2 


15969139200 


750 


29 


1 


E 32 3C 54 5A 66 


53 


16 


111122 4 


118540800 


20 


4 


3 


C 30 55 5A 66 69 


54 


24 


111122 6 


9676800 


5 


1 


7 19 2A 34 4C 52 61 


55 


32 


111222 4 


151200 


1 


1 


F 33 3C 55 5A 66 69 


Total: 


562949953421312 


33642660 


656103 
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Table A.3 
Continued 



^8 


The number of 






The SNF-class 




det 


SNF 


matrices 


7r-classes 


</}-ciasses 






representative 








0000000 





1 


1 


1 


n 
U 














1 





0000000 


1 


65025 


64 


16 














1 


2 





0000001 


1 


336954750 


2800 


295 














12 


3 





0000011 


1 


396683821800 


110064 


5758 














12 4 


4 





0000011 


2 


362237400 


596 


52 














3 5 6 


5 





0000111 


1 


144191294561160 


3696215 


114651 














12 4 8 


6 





0000111 


2 


801119894400 


65292 


2744 














1 6 A C 


7 





0000111 


3 


5797968120 


1422 


95 














3 5 9 E 


8 





0001111 


1 


17559952974446400 


88462953 


1874266 











1 


2 4 8 10 


9 





0001111 


2 


379659804531840 


3806686 


96326 











1 


2 C 14 18 


10 





0001111 


3 


11124304309440 


199374 


6235 











1 


6 A 12 1C 


11 





0001111 


4 


1070841280320 


27679 


957 











3 


5 9 11 IE 


12 





0001111 


5 


50409475200 


2430 


108 











3 


5 E 16 19 


13 





0001112 


2 


350465875200 


8104 


265 











3 


C 15 16 19 


14 





0011111 


1 


669716034190338240 


1150627540 


18733404 








1 


2 


4 8 10 20 


15 





0011111 


2 


46673270510307840 


114940091 


2053226 








1 


2 


4 18 28 30 


16 





0011111 


3 


3467174719659840 


11464276 


226089 








1 


2 


C 14 24 38 


17 





0011111 


4 


744008322193920 


3048569 


63679 








1 


6 


A 12 22 3C 


18 





0011111 


5 


90294382946880 


474025 


10919 








1 


6 


A 1C 2C 32 


19 





0011111 


6 


45671213145600 


258700 


5870 








3 


5 


9 16 2E 31 


20 





0011111 


7 


4508891956800 


33971 


883 








3 


5 


9 IE 2E 31 


21 





0011111 


8 


1853237232000 


15033 


409 








3 


5 


E 19 29 36 


22 





0011111 


9 


225909129600 


2474 


74 








3 


D 15 1A 26 39 


23 





0011112 


2 


244678270233600 


953097 


19406 








1 


6 18 2A2C 32 


24 





0011112 


4 


905427331200 


6401 


185 








3 


C 


15 1A 26 29 


25 





0011122 


2 


37302249600 


348 


18 








7 19 1E2A2D 33 


26 





0111111 


1 


5946448529329701120 


5204144555 


71348129 





1 


2 


4 


8 10 20 40 


27 





0111111 


2 


1255541169460515840 


1268895126 


18075782 





1 


2 


4 


8 30 50 60 


28 





0111111 


3 


201557577515938560 


230051193 


3411083 





1 


2 


4 18 28 48 70 


29 





0111111 


4 


78772224791393280 


98524334 


1498752 





1 


2 


C 


14 24 44 78 


30 





0111111 


5 


16904181842714880 


23305541 


367457 





1 


2 


C 14 38 58 64 


31 





0111111 


6 


15287255687531520 


21926711 


345484 





1 


6 


A 


12 2C 5C 62 


32 





0111111 


7 


2483072124099840 


3947501 


65221 





1 


6 


A 


12 3C 5C 62 


33 





0111111 


8 


1831786505418240 


3014954 


50240 





1 


6 


A1C 32 52 6C 


34 





0111111 


9 


601674419243520 


1055591 


17966 





1 


6 


1A 2A 34 4C 72 


35 





011111110 


451640549606400 


806664 


13584 





3 


5 


9 IE 2E 4E 71 


36 





011111111 


94403115878400 


185595 


ozyo 


n 
u 


3 


5 


9 IE 31 51 6E 


37 





011111112 


118893204864000 


230334 


4033 





3 


5 


E 


16 39 59 66 


38 





011111113 


19770573312000 


42957 


793 





3 


5 


E 19 36 56 69 


39 





011111114 


18696085632000 


39443 


695 





3 


5 


19 29 36 4E 71 


40 





011111115 


5859844300800 


12842 


235 





3 


D 


15 26 38 5E 61 


41 





011111116 


3094524518400 


7404 


157 





3 


C 


15 36 39 5A 65 


42 





011111117 


526727577600 


1380 


27 





3 


D 16 2E 39 5A 65 


43 





0111111 


L8 


702303436800 


1376 


24 





7 


19 2A 


34 4C 53 65 
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Table A.3 
Continued 



^8 



The number of 





det 


SNF 


matrices 


7r-classes 


0-classes 


representative 


44 





0111112 2 


25998918420787200 


31764328 


476566 





1 


2 C 30 54 58 64 


45 





0111112 4 


904945412044800 


1434616 


23722 





1 


6 18 2A 34 4C 52 


46 





0111112 6 


38868105830400 


70882 


1251 





3 


5 19 29 36 4E 51 


47 





0111112 8 


1777705574400 


4348 


114 





3 


C 31 55 5A 66 69 


48 





011111210 


175575859200 


408 


10 





7 


19 2A 34 4C 52 63 


49 





0111113 3 


70619902617600 


117219 


1972 





3 


5 18 28 49 4E 71 


50 





0111113 6 


181062604800 


385 


10 





3 ID 2D 36 3A 4E 71 


51 





0111114 4 


362125209600 


902 


23 





3 


C 35 3A 55 66 69 


52 





0111122 2 


40931905996800 


65455 


1203 





1 


E 32 3C 54 5A 66 


53 





0111122 4 


192036096000 


606 


26 





3 


C 30 55 5A 66 69 


54 





0111122 6 


12541132800 


66 


5 





7 19 2A 34 4C 52 61 


55 





0111222 4 


195955200 


5 


1 





F 33 3C 55 5A 66 69 


56 


1 


1111111 1 


3766962568171582080 


2363927011 


29610494 


1 


2 


4 8 10 20 40 80 


57 


2 


1111111 2 


3107221856321587200 


1958051993 


24598561 


1 


2 


4 8 10 60 AO CO 


58 


3 


1111111 3 


1128344550375409920 


717105693 


9069162 


1 


2 


4 8 30 50 90 E0 


59 


4 


1111111 4 


798113338051276800 


508274669 


6438161 


1 


2 


4 18 28 48 88 F0 


60 


5 


1111111 5 


280558398045864960 


180518667 


2303868 


1 


2 


4 18 28 70 B0 C8 


61 


6 


1111111 6 


391839981330309120 


249971410 


3172566 


1 


2 


C 14 24 58 B8 C4 


62 


7 


1111111 7 


92717618729258880 


60152437 


772896 


1 


2 


C 14 24 78 B8 C4 


63 


8 


1111111 8 


98081405804067840 


63280071 


811462 


1 


2 


C 14 38 64 A4 D8 


64 


9 


1111111 9 


46392962843324160 


30131275 


388116 


1 


2 


C 34 54 68 98 E4 


65 


10 


111111110 


49370839378882560 


31886563 


408631 


1 


6 


A 12 3C5C 9C E2 


66 


11 


111111111 


14214384381012480 


9334363 


121275 


1 


6 


A 12 3C 62 A2DC 


67 


12 


111111112 


25287474431600640 


16375045 


211035 


1 


6 


A1C2C 72 B2CC 


68 


13 


111111113 


6076097931697920 


4012940 


52419 


1 


6 


A 1C 32 6CAC D2 


69 


14 


1 1 1 1 1 1 1 14 


8661618203857920 


5648182 


72918 


1 


6 


A 32 52 6C 9C E2 


70 


15 


111111115 


4660876422921600 


3060154 


39958 


1 


61A2A4C 70 BC C2 


71 


16 


111111116 


3363977985177600 


2214742 


29108 


1 


6 


18 2A 6C 72 B4 CA 


72 


17 


111111117 


1206477746611200 


804732 


10608 


1 


6 1A 2C 5C 72 B4 CA 


73 


18 


111111118 


2308713728025600 


1503637 


19409 


1 


E 


32 54 68 98 A6 CA 


74 


19 


111111119 


548565282316800 


367458 


4865 


3 


5 


9 31 51 6E 9E El 


75 


20 


1 1 1 1 1 1 1 20 


883792178841600 


583710 


7689 


3 


5 


E 19 36 68 A9 D6 


76 


21 


111111121 


420654153830400 


280261 


3719 


3 


5 


E 19 36 69 A9 D6 


77 


22 


1 1 1 1 1 1 1 22 


358862424883200 


237040 


3085 


3 


5 


E 32 56 69 99 E6 


78 


23 


1 1 1 1 1 1 1 23 


113325986995200 


76671 


1023 


3 


5 


E 36 56 69 99 E6 


79 


24 


1 1 1 1 1 1 1 24 


228366581990400 


152796 


2117 


3 


5 


18 68A9B6CE Dl 


80 


25 


1 1 1 1 1 1 1 25 


59747916211200 


40146 


532 


3 


5 


18 2E 69 76 Bl CE 


81 


26 


1 1 1 1 1 1 1 26 


73204159795200 


48393 


630 


3 


5 


19 2E 69 76 Bl CE 


82 


27 


1 1 1 1 1 1 1 27 


33155489203200 


22046 


293 


3 


5 


39 59 6E 76 9E El 


83 


28 


1 1 1 1 1 1 1 28 


34709152665600 


23559 


324 


3 


C 


15 36 5A 65 B9 C6 


84 


29 


1 1 1 1 1 1 1 29 


9118971187200 


6180 


82 


3 


C 


31 55 7A 96 D9 E5 


85 


30 


1 1 1 1 1 1 1 30 


28015323033600 


18147 


238 


3 


c 


31 54 6A 9A A6 C9 


86 


31 


111111131 


3621252096000 


2436 


32 


3 


D 


16 2E 5A 75 B9 C6 



The SNF-class 
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Table A.3 
Continued 



-4s 


The number of 




Ine bNt -class 
representative 


det 


SNF 


matrices 


7r-classes 


^-classes 


87 


32 


111111132 


5423648025600 


3953 


68 


3 C 31 56 6A 9A A6 C5 


88 


33 


1 1 1 1 1 1 1 33 


2806470374400 


1866 


27 


3 D IE 35 66 79 AA D5 


89 


34 


1 1 1 1 1 1 1 34 


2831160729600 


1815 


23 


3 D 31 54 6A 9A A6 C9 


90 


35 


111111135 


757170892800 


491 


7 


3 ID 2E 56 69 9A Bl CD 


91 


36 


1 1 1 1 1 1 1 36 


1327792435200 


873 


13 


3 D 35 56 69 99 AE C5 


92 


37 


111111137 


131681894400 


81 


1 


7 19 2A4C 71 A5 CB D6 


93 


38 


1 1 1 1 1 1 1 38 


592568524800 


389 


5 


3 ID 2E 56 69 9A B5 CD 


94 


39 


1 1 1 1 1 1 1 39 


65840947200 


45 


1 


7 19 2A 56 6D 9C B3 CB 


95 


40 


1 1 1 1 1 1 1 40 


263363788800 


200 


4 


7 19 2A 4C 71 96 AD CB 


96 


42 


1 1 1 1 1 1 1 42 


65840947200 


45 


1 


7 19 3E 63 AA B5 CC D2 


97 


4 


1111112 2 


264489939127895040 


166494351 


2095861 


1 2 4 18 60 A8 BO C8 


98 


8 


1111112 4 


48628694487582720 


30927272 


394107 


1 2 C 30 54 68 98 A4 


99 


12 


1111112 6 


8317470133324800 


5286249 


67857 


1 6 A 32 52 6C 9C A2 


100 


16 


1111112 8 


1671575454259200 


1087782 


14500 


1 6 18 62 AA B4CC D2 


101 


20 


111111210 


300128743756800 


194974 


2641 


1 E 32 54 68 98 A4 C6 


102 


24 


111111212 


119718045619200 


80651 


1172 


3 5 18 60 A9 B6 CE Dl 


103 


28 


111111214 


14388990336000 


9935 


157 


3 C 15 3A 65 A5 D6 D9 


104 


32 


111111216 


4180900147200 


3259 


65 


3 C 31 56 6A 9A A6 CI 


105 


36 


111111218 


1360712908800 


1103 


23 


3 D 31 54 6A 9A A6 CI 


106 


40 


1 1 1 1 1 1 2 20 


271593907200 


232 


5 


3 D 35 59 6E 9E A9 C5 


107 


44 


111111222 


76814438400 


74 


2 


3 ID 2E 56 79 9A B5 CD 


108 


48 


1 1 1 1 1 1 2 24 


16460236800 


28 


1 


7 19 3E 61 AB B5 CC D2 


109 


9 


1111113 3 


5728974559056000 


3662516 


47056 


1 6 A 30 50 92 9C E2 


110 


18 


1111113 6 


292630089830400 


187420 


2455 


1 6 3A5A 6C 74 9C E2 


111 


27 


1111113 9 


6466129689600 


4234 


71 


3 C 31 54 7A 9A A6 C9 


112 


36 


111111312 


285310771200 


219 


7 


3 1C 65 6A A6 Bl C9 D2 


113 


45 


111111315 


2743372800 


6 


1 


7 39 5A 6C 9C AB B6 Dl 


114 


16 


1111114 4 


282699080294400 


186028 


2556 


1 6 18 6A 74 AA CC D2 


115 


32 


1111114 8 


724250419200 


713 


21 


3 C 35 3A 56 69 99 A6 


116 


25 


1111115 5 


2339411155200 


1497 


22 


3 D 15 26 5E 61 B8 C5 


117 


36 


1111116 6 


142655385600 


136 


4 


7 19 2A4B 74 8C D2 El 


118 


8 


1111122 2 


2260349894476800 


1421783 


18397 


1 2 1C 64 78 A8 B4 CC 


119 


16 


1111122 4 


136245037824000 


90153 


1346 


1 6 18 60 AA B4 CC D2 


120 


24 


1111122 6 


6530011084800 


5175 


114 


1 E 32 54 68 98 A4 C2 


121 


32 


1111122 8 


R9^/issQQS/inn 
ozo4ooyyo4uu 




15 


3 D 31 55 oA 9A Ab CI 


122 


40 


111112210 


54867456000 


96 


5 


3 D 31 56 6A 9A A6 CI 


123 


48 


111112212 


16460236800 


49 


3 


3 1C 65 7A A9 B6 CE Dl 


124 


56 


111112214 


391910400 


6 


1 


3 ID 65 7A A9 B6 CE Dl 


125 


32 


1111124 4 


98761420800 


134 


4 


3 1C 64 78 A9 B2 CA Dl 


126 


27 


1111133 3 


101504793600 


96 


6 


7 19 2A 4B 74 8D Bl D6 


127 


16 


1111222 2 


1082717798400 


754 


17 


31C 64 78 A9AA B5CD 


128 


32 


1111222 4 


28217548800 


74 


4 


1 IE 66 78 AA B4 CC D2 


Total: 


18446744073709551616 


14685630688 


199727714 
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Table A.4 

Transposed incidence matrix M-j containing all M n , 1 < n < 7. Symbol at position 
(s', s) carries information about m SjS > (7): •, *, o denotes respectively 1, explained 
by Lemma 10, and not explained by Lemma 10. 



000 
000 

111 
111 
111 



uuouu 

00000 

11111 
11111 
11111 
11111 



000000000 

111111111 
111111111 
111111111 
111111111 
111111111 



uuu 
111 
111 
111 

112 

222 



111111111111111111 
111111111111111111 
111111111111111111 
111111111111111111 
111111111111111111 
111111111111111111 



TTTTT 
11111 
11111 
11111 
11111 
22222 



000 
123 



00000 
12345 



000000000 
123456789 



000 

242 



000000000111111111 
123456789012345678 



00001 
24680 



0000000 



0000000 1 



0000001 1 



0000011 1 
0000011 2 



0000111 1 
0000111 2 
0000111 3 



••• 

*•* 
**• 



0001111 
0001111 
0001111 
0001111 
0001111 



0001112 2 



••• 
••• 
••• 
••• 



*•*•* 

**•** 
***** 
***** 



*•* 



***** 



0011111 
0011111 
0011111 
0011111 
0011111 
0011111 6 
0011111 7 
0011111 8 
0011111 9 



0011112 2 
0011112 4 
0011122 2 



•** 
••* 
•*• 
••* 

•** 
• oo 
o** 
••* 
•*o 



• •••• 

• •••• 

• •••• 



••••••••• 

*•*•*•*•* 

**•**•**• 
***•***•* 

********* 
********* 
********* 
********* 
********* 



• *o 

• •• 

*** 
o*o 

*** 
*** 
*** 
*oo 
*** 



••* 
••* 

*•* 



*•*•* 
*•*•* 

***** 



********* 
********* 
********* 



• •• 

*»o 
**• 



0111111 1 
0111111 2 
0111111 3 
0111111 4 
0111111 5 
0111111 6 
0111111 7 
0111111 8 
0111111 9 
0111111 10 
0111111 11 
0111111 12 
0111111 13 
0111111 14 
0111111 15 
0111111 16 
0111111 17 
0111111 18 



0111112 2 
0111112 4 
0111112 6 
0111112 8 
0111112 10 



***** 
••*** 
•*•** 
••*•* 

•***• 
•••** 

***** 
••*•* 

•*•** 

•o**o 
o**** 

• ••o* 
o**** 
•o*** 
•*o*o 

• •*o* 
o**** 

• •o** 



••••••••• 

••••••••• 

••••••••• 

••••••••• 

••••••••• 



••••••••• 



*** 
••* 
*** 

• •o 
*** 
••* 
*** 

• •o 
*** 
••* 



*** 



• •o 

*** 
••* 
*** 

• •o 



*** 



*•*•*•*•*•*•*•*•*• 

**•**•**•**•**•**• 
***•***•***•***•** 
****************** 
****************** 



****************** 
****************** 
****************** 
****************** 
****************** 
****************** 
****************** 
****************** 
****************** 
****************** 
****************** 



• •••• 

***** 
o*o*o 

***** 
***** 

***** 
*o*»* 
***** 
***** 



***** 
**o** 
***** 
***** 
***** 
***o* 
***** 
***** 



0111113 3 

0111113 6 

0111114 4 
0111122 2 
0111122 4 
0111122 6 
0111222 4 



• **o* 
***** 

• **o* 

• **o* 

oo*oo 



*•*•*•*•* 
*•*•*•*•* 
*•*•*•*•* 
*•*•*•*•* 
*•*•*•*•* 



• •• 

• •• 

• •• 
••• 
••• 



****************** 



****************** 
****************** 



****************** 



• •••• 

*•*•* 

***** 
***** 
***** 



•*o** 

• •o** 

• •*o* 
*»*o* 
*»*o* 
*o*o* 

***** 



**•**•**• 

**•**•**• 

***•***•* 

********* 
********* 
********* 
********* 



*** 
*** 

o»o 



••• 

**• 



****************** 

****************** 
****************** 
****************** 
****************** 
****************** 

****************** 



***** 
***** 
***** 
***** 
***** 
***** 
***** 
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Table A.4 



Continued. 



s 

a' 

























1 


01) 

00 
00 
00 

11 
11 


000 

000 
000 

111 
111 
111 


00000 
00000 

11111 
11111 
11111 
11111 




1 
1 
1 
2 


000000000 
111111111 
111111111 
111111111 
111111111 
111111111 


000 

111 
111 
111 

112 

222 


111111111111111111 
111111111111111111 
111111111111111111 
111111111111111111 
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The part of 5g corresponding to nonsingular part of Ag. 
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Table A.6 

The part of the transposed incidence matrix M§, corresponding to the regular part 
of Sg. Symbol at position (s',s) carries information about m S)S i (7): •, *, o denotes 
respectively 1, explained by Lemma 10, and not explained by Lemma 10. 
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Table A.6 
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Table A.7 



Lower bounds for a n and matrices with high extension spectra, 10 < n < 19. 



n 


a n > 




det A n -i 






10 


259 


2 


110 


257 


[7, 39, 5A, 9C, 51, 149, 174, 193, 1AA] 


11 


739 


6 


291 


679 


[F, 71, 56, 13A, 1C3, IDC, 256, 299, 2EC, 325] 


12 


2107 


19 


779 


1894 


[F, 73, 195, IE A, 2A6, 35C, 456, 535, 565, 
659, 703] 


13 


6157 


18 


2201 


5618 


[IF, 53, 17C, 3A5, 649, 656, 732, A6E, AB8, 
553, C35, D8E] 


14 


19073 


40 


6731 


16821 


[3F, 1C7, 2D9, 76A, CAD, CF2, 594, 1575, 
1685, 195,4, 19A9, 1AU, 1E13] 


15 


58741 


46 


23288 


53117 


[7D, 38F, 552, 555, 1895, 1954, 1529, 2AEA, 
2D1C, 3254, 3353, 34C9, 3C27, 1645] 


16 


185693 


190 


67832 


161599 


[FD, 715, 553, 1529, 3245, 3652, 3995, 5370, 
55C6, 5,49.4, 61AB, 6C53, 6524, 27C8, 2975] 


17 


610187 


480 


213175 


517794 


[155, 535, 17C6, 3,453, 6495, 6564, 7325, 
,4652, AB8D, 5535, C356, D8A7, 5C49, 4591, 
725C, 599,4] 


18 


2039033 


697 


709503 


1719277 


[359, 1C75, 2595, 76AC, C452, C527, 5949, 
15755, 16855, 195,43, 19A9C, 1,4645, 15135, 
1458,4, 17,433, 33C6, 1,4570] 


19 


6478579 


54 


2331887 


4663774 


[759, 3857, 5513, 5955, 19277, 15599, 1CC45, 
29585, 2,4531, 254C5, 30556, 36295, 37125, 
1354C, 1EBC2, 35858, 2546,4, 5754] 
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